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ABSTRACT 



Aims. We have studied the characteristic timescales of 80 AGNs at 22, 37 and 90 GHz examining the properties of the wavelet method 
and comparing them to traditional Fourier-based methods commonly used in astronomy. 

Methods. We used the continuous wavelet transform with the Morlet wavelet to study the characteristic timescales. We also gain 
information when the timescale is present in the flux curve and if it is persistent or not. 

Results. Our results show that the sources are not periodic and changes in the timescales over a long time are common. The property 
of wavelets to be able to distinguish when the timescale has been present is superior to the Fourier-based methods. Therefore we 
consider it appropriate to use wavelets when the quasi-periodicities in AGNs are studied. 

Key words, galaxies: active - methods: statistical - radio continuum: galaxies 



1. Introduction 

Wavelet methods have been used for the past 10 years to study 
the variab ility timescales in different types of astron omical ob- 
jects (e.g. lFosteij[T996l: iPriestlevll 1 997t IScargl^l 1 997 ). Wavelets 
have properties superior to Fourier based methods, the most use- 
ful being the locality of the timescale. In Fourier based meth- 
ods a single sinusoid is fitted to the whole light curve, whereas 
wavelets enable us to study temporally local timescales or pe- 
riods that are present in the light curves. This is a very good 
property when studying flux curves of AGNs, which are almost 
never strictly periodic over long times. 

Wavelets were used to study vari ability timescales of a 
larger sample of AGNs by iKellv et al.l ((2003). They used the 
cross-wavelet transform to study quasi-periodic variations in the 
Pearson-Readhead VLBI sample, monitored by the University of 
Michigan Radio Observatory (UMRAO). A quasi-periodic be- 
haviour was found in a little over half of the 30 sources they 
studied. Individual sources have also been studied f or quasi- 
periodicity and periodicity using wavelets (e.g. iHughes et al.l 
ll998HLehtollI999tlKadIer et al.ll2006tlciprini et alJl2007h . 

In iHovatta et al.l (l2007h . hereafter Paper I, we studied the 
variability timescales of a large sample of AGNs using the struc- 
ture function (SF), the discrete autocoiTelation function (DCF) 
and the Lomb-Scargle periodogram. We found that larger flares 
in these sources happen quite rarely, only once in 4-6 years 
while short-term variations are continuously present in the flux 
curves. Also, many of the sources had changed their behaviour 
during the 25 years of monitoring which is one reason why we 
wanted to study the sample again using wavelets. In Paper I we 
also studied the properties of the methods and their differences 
by searching for the characteristic timescales. We noticed that 
there are several problems when using periodograms to study 
the timescales because none of the sources in our sample showed 



strict periodicities in their flux curves or were well characterised 
by sinusoids. 

In this paper we will compare the results from this wavelet 
analysis to Paper I. This paper is organised as follows: In Sect.|2] 
we describe the source sample and data used. Section[3]includes 
the description of the method and in Sect. |4] we present the re- 
sults. The discussion is in Sect.|5]and the conclusions are drawn 
in Sect. |6] 

2. The sample and observations 

We used the same sample of 80 AGNs selected from the 
Metsahovi monitoring list that was used in Paper I. The sample 
consists of 24 BL Lacertae Objects (BLOs), 23 Highly Polarised 
Quasars (HPQs), 28 Low Polarisation Quasars (LPQs) and 5 
Radio Galaxies (GALs). The five sources with no information 
about their optical polarisation are considered to be LPQs in this 
study and are included in the numbers above. 

All sources have been monitored for at least 10 years at 22 
or 37 GHz. These are all bright sources with a flux density of 
at least 1 Jy in the active state in at least one of the two fre- 
quencies. We used the 22, 37 and 87 GHz data obtained with 
the Met sahovi 14-metre radio telescope over 25 years of mon- 



itoring (*SaIo nen et al 



2005; Nieppo la et aP 



19871: iTerasranta et alJ[T992ll 19981 12004 
2007h . The 37 GHz data of HPQs, LPQs 



and GALs from December 2001 until March 2005 are partially 
unpublished. The details of the observation method and data re- 
duction processes are given in ITerasranta et al.! (Il998b . In ad- 
dition we used 90 GHz data obtained with the Swedish-ESO 
Submillimetre Telescope (SEST) at La Silla, Chile, which was 
used in our monitoring campaign to sample the high frequency, 
90 and 230 GHz variability of southern and equatorial sources 
(iTornikoski et al.l 11996). The monitoring campaign at SEST 
lasted from 1987 to 2003 from which the data from 1994.5-2003 
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are unpublished. Complementary high frequency data at 90 GHz 
were also collected from the literature ("Steppe et £111198811 19921 
[1993; Reuter et al. 1997). The median intervals between the ob- 
servations at 22, 37 and 90 GHz are 31, 41 and 47 days, respec- 
tively. The values depend on the source and at 37 GHz the min- 
imum average sampling rate was 6.8 days for the source 3C 84 
which is used as a secondary calibrator in the Metsahovi obser- 
vations. The maximum average value at 37 GHz was 186.4 days 
for the source 2234-1-282. 



3. The wavelet method 

We have analysed the light curves using Morlet wavelets. They 
can be understood as local wave packets. The method is closely 
related to chirp analysis methods. Short Time Fourier Transform 
methods, Gabor transforms and also to the sonogram or spectro- 
gram analysis used widely in audio signal processing. In wavelet 
analysis the length of the analysis "package", the kernel, is pro- 
portional to the timescale of interest, or in terms of frequency, 
proportional to the inverse of the frequency. We have used a 
Morlet wavelet, which can be understood mathematically as a 
localised plane wave tapered by a Gaussian function. For longer 
timescales of interest we thus use lower frequency sinusoids 
than for higher frequencies according to the following equations. 
Wavelet transforms are functions of both frequency, /, (or the re- 
spective timescale) and the local time, r. 

g'if, T) = exp ^-icfit - T) - i(/(f - T))2j , (1) 

where we adopt c = 2n, which characterises the amount of ta- 
pering. The power of the wavelet transform, known as the scalo- 
gram, is defined for evenly spaced data as 

W{f,T)^f-{sHf,T) + CHf,T)), (2) 

where 

S (/, t)^Yj (2^/(0 - ^)) exp I — (/(f; - T))' 1 , (3) 



C{f, r) = J mj cos [infitj - r)) exp |^(/(f,. - t)^^ , (4) 

where nij and tj are respectively the flux density and the time of 
the observation j. 

For unevenly spaced data we have to weight each point with 
the inverse of the local density of points. The locality, a property 
of wavelets and other similar methods mentioned above, is in 
strong contrast to global methods such as Fourier transforms or 
structure functions, which treat the whole light curve as one en- 
tity. A change in a few critical points can change the final result 
in global methods rather significantly, whereas in wavelet anal- 
ysis and other local methods a chance point in the data will only 
affect the analysis locally. As we gain information about the lo- 
cal nature of the signal we lose in the resolution with which we 
can determine the timescale (or frequency) of the signal. The 
length of the kernel used limits our resolution to about 0.07 in 
log space, which is significantly worse than what would be ob- 
tained e.g. from the Fourier analysis of a long lasting sinusoid 
with a constant period. Wavelet transforms with longer kernels 



could also be applied to the data as intermediate forms approach- 
ing Fourier transforms in accuracy, but in doing so one loses the 
locality information according to 

AfAv > 6, (5) 

which is analogous to the Heisenberg uncertainty principle. Here 
6 can be considered as unity if we are interested in separating 
two nearby peaks. If we rather wish to measure the location of 
a single signal in weakly noisy data then 6 is of the order of 
yfiTJWj, where is the local number of points. 

First we calculate the scalogram for our unevenly spaced 
data. We then calculate a scalogram of normally distributed 
Gaussian white noise with a variance of unity sampled as was the 
original data. By repeating this 20 times we get a good estimate 
of the scalogram of this white noise. We then normalise our orig- 
inal scalogram by the white noise scalogram. If our data were 
pure noise, this noise-normalised scalogram would have every- 
where an expected value equal to the variance of our measure- 
ments. Alternatively, we argue that by taking the square root of 
the noise-normalised scalogram and dividing it by the (known) 
measurement errors we get an estimate of the significance of the 
variations detected. If a variability timescale has a variance of 
> 9 as expected from measurement errors, we consider that a 
significant detection and if the variance is between 4 and 9 times 
the variance of measurements, we consider that a marginal case 
(marked as a weak timescale in Table [T]i. These numbers were 
selected to represent "2 and 3 sigma" levels, but we caution that 
these should be considered indicative only, as we have not yet 
carried out full simulations of the probability density distribu- 
tions of the noise properties in noise-normalised scalograms. 

4. Wavelet timescales 

4.1. Interpretation 

The advantage of wavelets is that in addition to finding the vari- 
ability timescale in the frequency domain, we can see at what 
moments it has been present. This enables us to see whether the 
timescale is present during the whole observed time or is tran- 
sient. In addition we can identify typical flare timescales i.e. flare 
duration or rise and decay times and notice whether one large 
outburst dominates the timescale obtained. Also short-term vari- 
ability is easily detected. 

We searched for both long-term variability, which is present 
most of the time in the flux curves, and typical timescales in 
which the flares occur. For the long-term timescales we first 
identified the timescale from the frequency axis and then deter- 
mined how long it has been present on the time axis. This way 
we can study how many times the timescale has repeated during 
the time it has been present and also how long time it has been 
present compared to the total observing period. Figure [T] shows 
the wavelet transform for the source 1156+295 (4C 29.45). We 
find a timescale of 3.4 years (lO^'^^* years) which is present in 
the latter half of the flux curve, starting at 1995. The timescale 
is marked with a circle in the plot. Furthermore we identified 
typical flare timescales in the cases where more than one flare 
occurs on approximately the same timescale. For 1 156+295 it is 
seen at 1.7 years (lO"'^-' years). This timescale is also shown in 
the plot as circles at the beginning of the time series at the time 
of the first two flares and later simultaneously with the long-term 
trend. The timescale can be identified with the durations of these 
flares which are of about 1-2 years. 

Often the long-term timescales are either slowly shortening 
or increasing in time which is seen as rising or declining trends in 
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1985 1990 1995 2000 



Fig.l. Wavelet transform of the source 1156+295 at 22 GHz. 
The y-axis is frequency 1/year so that a timescale of 1 year 
corresponds to frequency and a timescale of 10 years corre- 
sponds to frequency -1 in the log-axis. A long-term variability 
timescale of 3.4 years (lO"'^-' years) is visible in the latter half 
of the flux curve starting at 1995. A flare timescale of 1.7 years 
(^lQ0-2i years) is also seen. 



the wavelet plot. These are often connected to flares and usually 
the timescales shorten towards the peak of a large flare. 

4.2. Long-term trends 

We were able to determine a long-term timescale for 122 cases 
out of the total 192 cases for which we calculated the wavelet 
transform. In addition to the timescale, we determined how long 
it is present in the flux curve. This way we were able to calculate 
how many times the cycle has repeated itself. The timescales 
are listed in Table [1] where the B1950-name, a well-known 
alias (when applicable), classification, total monitoring time, 
timescale and number of cycles and timescale of flares are tabu- 
lated for each frequency band. If there is more than one timescale 
present the most significant one is given first. In 18 sources (~l/4 
of the sample) we found long-lasting cycles, continuing for at 
least four cycles. Five of these sources showed cycles at two 
separate frequencies and three sources at all the three frequency 
bands. There were 31 cases where the variability timescale did 
not last for two cycles. In Table |2] we have calculated averages 
and medians for the timescales continuing for at least two cycles. 
They are given for each frequency and source class separately. 
In addition we corrected the timescales for redshift and these 
averages are also shown in Table |2] 

We ran a Kruskal-Wallis analysis to study whether the 
timescales at different frequencies and source classes differ from 
each other The timescales at 90 GHz differ from lower fre- 
quencies with a 95% confidence limit with the timescales be- 
ing shorter. The result can be affected by the smaller number 
of sources at 90 GHz and shorter monitoring period which also 
shortens the timescales obtained at 90 GHz. We studied the dif- 
ferences between the source classes at 22 and 37 GHz and at 
both frequencies the classes do not differ significantly from each 
other if the observational uncorrected timescales are considered. 
When we study the redshift corrected timescales we see signif- 
icant differences between the BLOs and LPQs at 37 GHz and 
in addition there is strong indication that also the HPQs and 
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BLOs differ. Similar results were obtained in Paper I, when pe- 
riodogram timescales were studied. 

We were also interested in sources in which the timescale is 
present during most of the flux curve. By visually extracting the 
duration of the timescale from the wavelet plot, we calculated 
how much of the total monitoring time it has been present. In 
27 cases the timescale was present at least 90% of the total ob- 
serving time. Seven of these, however, were cases in which the 
timescale had not repeated twice during the period it was visible. 

4.3. Flare timescales 

In addition to long-term trends we determined a flare timescale 
which describes the duration or rise and decay times of flares. 
These timescales are seen in the wavelet transform at the times of 
the outbursts. We determined such a timescale in 53 cases. They 
are also tabulated in Table [T] We have taken into account cases 
for which a similar timescale seems to be present during more 
than one flare and therefore represented a typical flare timescale 
for the source. The values we obtained vary between 0.34-2.17 
years at 22 and 37 GHz. The average value at 22 GHz is 1.16 
years and at 37 GHz it is 1 .03 years. At 90 GHz, there were only 
6 sources which showed repeating flares and the timescales were 
shorter, varying between 0.34 and 0.87 years, with an average of 
0.69 years. 

4.4. Individual sources 

Even though none of the sources studied here seems to show a 
strict periodicity in radio frequencies, many of them show quasi- 
periodic behaviour. We are interested in sources that have a vari- 
ability timescale lasting four cycles in at least two of the fre- 
quency bands. There are eight such sources and their properties 
are discussed in more detail below. The wavelet plots for these 
sources at 37 GHz are shown in Figs.|2]-|9](only available in elec- 
tronic form via http://www.aanda.org). 

DA 55 (0133+476): This HPQ source has been monitored in 
Metsahovi for over 20 years at 22 and 37 GHz and for almost 
15 years at 90 GH z. It has shown hi gh apparent pattern velocity 
of 18c in the jet jPiner et al.l l2007h and has been also a target 
of many other studies usin g Very Long Baseline Interferometry 
(VLBI) observations (e.g. ' ListeJl200iI iKeUermann et al.ll2004l: 
Lister & Homan 2005). It is also included in the wavelet study 
of IKellv et al.1 ( 120031) which uses the UMRAO 4.8, 8, 14.5 GHz 
monitoring data. They found a timescale of 2.3 years for the 
source at 14.5 GHz. This is very close to the timescale of 2.2 
years which we obtained for the source at both 22 and 37 GHz. 
At 90 GHz we only found a flare timescale of 0.7 years for 
the source. The long-term timescale in our analysis has a ris- 
ing trend so that it is shortening towards a peak of a large flare. 
The timescale is present only for about 40% of the flux curve 
at both 22 and 37 GHz but still has repeated 4.4 times at both 
of the frequency bands. In Paper I we could not find any DCF 
or periodogram timescale for the source. When the flux curve is 
examined we can clearly see that the timescale obtained here is 
related to a period when the source had many similar flares and 
the total flux level rises. This also explains the rising trend in the 
wavelet plot which is caused by more frequent individual flares 
towards the peak of the larger event. 

3C 120 (0430+052): This radio galaxy has been monitored 
in Metsahovi for over 20 years at 22 and 37 GHz and for 15 
years at 90 GHz. It is a well-studied object at all wavelengths and 
has been a target of many multiwavelength campaigns. Its ra- 
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Table 2. Averages and medians of long term timescales that have lasted for at least two cycles. Values are given for each frequency 
and source class separately. Also the redshift corrected averages are shown. 
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dio jet has been studied and modelled using VLBI observations 
by several authors (e.g. IWalker et al.ll200lb iGomez et al.ll200lt 
iHardee et al.ll2005"l) . Our wavelet analyses reveal a timescale of 

4.3 years at both 22 and 37 GHz. At 22 GHz the timescale has 
been present for 85% of the total observing period lasting for 

4.4 cycles. At 37 GHz it has been present for 84% of the time 
and persisted for 4.6 cycles. In Paper I the DCF analyses showed 
a timescale of 4.2 years at 22 GHz which is very close to what 
we obtained here. Both frequencies also show flare timescales 
of 1.4 years and in addition 0.5 years at 22 GHz and 0.3 years at 
37 GHz. These timescales are also seen in the DCF and SF anal- 
yses of Paper I. At 90 GHz the timescale obtained is 2.7 years 
and it has been present for 73% of the time and continued for 4 
cycles. This timescale is also very close to the DCF timescale of 
2.9 years. 

Visual inspection of the flux curve at 37 GHz shows that in- 
deed many of the larger flares have approximately 4.5 to 5 years 
between them but there are also additional flares in between. The 
last big flare included in our analysis at 37 GHz peaked in the au- 
tumn of 2003. The next big flare was observed at the end of 2006 
with an interval of 3 years between the flares. This indicates that 
the source is not strictly periodic but has shown a characteristic 
timescale of 4.3 years. 

0736+017: This HPQ source was not detected by EGRET 
but has been listed as a possible y-ray source to be detected 
with t he new satellite missions AGILE and GLAST (iBach et al.l 
I2007h . In Metsahovi it has been monitored for over 20 years 
at 22 and 37 GHz and for 1 1 years at 90 GHz. It exhibits rapid 
variability and i n the optical dom ain it has been stud ied by sev- 
eral authors (e.g.l Clements et alT2003 ; Ramfrez et al.t2004 ). The 
wavelet plot has a complex structure at all the frequency bands 
and we found multiple timescales for the source. The most sig- 
nificant and long-lasting ones are 2.7, 3.4 and 1.7 years at 22, 37 
and 90 GHz, respectively. The timescales have repeated them- 
selves 4 times at all the frequency bands and are also present for 
about half of the flux curve. The DCF analysis of Paper I shows 
very similar timescales of 2.8 years at 22 and 37 GHz and 1.2 
years at 90 GHz. At 22 and 37 GHz there are multiple timescales 
present also in the DCF which confirms the complex variability 
of this source. Therefore it is also difficult to predict when the 
next flares could be expected. 

OJ 248 (0827+243): This radio quasar was found by 
EGRET to be a bright y-ray source. It has been monitored at 
Metsahovi for about 11 years both at 22 and 37 GHz. It has 
also b een a target of multiwavelength campaigns (e.g. lBach et alj 
l2007h and has shown extremely fast apparent pattern velocity 
exceeding 25c in the jet dPiner et al. 2006) . The kiloparsec- scale 
jet is studied in more detail in Tjorstad & Marscheii (l200l . We 



found a timescale of 1 .4 years at both frequency bands but it ap- 
pears quite weak in strength and is classified as a weak timescale 
at 22 GHz. At 22 GHz it has been present for 54% of the time 
and is discontinuous. It repeated 4.2 times during the time it has 
been present. At 37 GHz the timescale has been present for 77% 
of the time and has a shortening trend changing the timescale 
from 1 .4 to 1.1 years. The timescale repeated 6.3 times. In Paper 
I we measured timescales of 1 .7 and 1 .4 years at 22 and 37 GHz 
respectively from the DCF analyses. 

OJ 287 (0851+202): This BLO type object is one of the 
most studied AGNs. In Metsahovi it has been monitored for over 
25 years at 22 and 37 GHz and at 90 GHz for 15 years. In the 
optical there is strong evidence of a period of about 12 years 
which is explai ned b y several binary black ho l e models (e.g. 
Sillanpaa et a ll 119881: iLehto & Valtoneiil 1 1 9961: IValtaoia et al] 
200Q; iValtonen et al.l 12008*). We can neither confirm nor dis- 
prove the existence of the same period at our frequency bands 
because the monitoring time is too short compared to the optical 
historical light curve which is over 100 years l ong. The sourc e 
has also been studied with wavelets by Hugh es et al.l ( Il998l) . 
They used the continuous wavelet transform to study more than 
20 years of data obtained in the University of Michigan Radio 
Observatory (UMRAO) at 4.8, 8 and 14.5 GHz. They found a 
persistent timescale of 1 .66 years in the total flux and polarisa- 
tion and another timescale of 1.12 years dominating the activity 
in the 1980s. They explained the modulations with a shock-in-jet 
model in which the permanent timescale is associated with the 
quiescent jet and the shorter timescale is due to the passage of a 
shock. 

We found a timescale of 1 .4 years at 22 and 90 GHz and 1 .7 
years at 37 GHz. Both of thes e can be considere d to be very sim- 
ilar to the ones obtained by Hugh es et alj (Il998l) . given the ambi- 
guity in determining the exact time scale in our analysis. A very 
interesting phenomenon seen in the wavelet transform is that the 
timescale either weakened or totally disappeared between 1993 
and 2000. Due to this the timescale is present only for 2/3 of 
the time in all the frequency bands. It can be clearly seen in the 
flux density curve in Fig.|6]that the flux density level is generally 
lower during that time, even though smaller flares occur contin- 
uosly. It is also possible to see in the wavelet plot that the num- 
ber of rapid flares is lower during that period than in the more 
active states. Nevertheless the timescale repeated 10.2 times at 
22 GHz, 9.6 times at 37 GHz and 7.2 times at 90 GHz. AH the 
timescales have a corresponding one in the DCF analysis of 
Paper I within 6 months of the one obtained from wavelet anal- 
ysis. When studying the flare timescales we noticed timescales 
of 0.4 and 0.3 years at 22 and 37 GHz, respectively. These agree 
with timescales from the SF analysis in Paper I. It is difficult to 
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associate the timescales with individual flares because the source 
exhibits continuous variabihty and is almost never in a quiescent 
state. 

OL 093 (1055+018) : This HPQ source has been observed in 
Metsahovi for over 20 years at 22 and 37 GHz and for almost 
13 years at 90 GHz. The j et of this source is suggested to have 
a spin e-sheath structure ( Attrid ge et al.l 119991: iPushkarev et al.l 
l2005h . We obtained a timescale of 4.3 years for the source at 
both of the 22 and 37 GHz frequencies. The timescale repeated 
4.4 times and is present for 85 and 80% of the flux curve at 22 
and 37 GHz, respectively. Very similar timescales are measured 
in Paper I, where a DCF timescale of 4.7 and 4.4 years are ob- 
tained at 22 and 37 GHz, respectively. At 90 GHz we obtain a 
timescale of 2.7 years but it repeated only 2.6 times. When ex- 
amining the flux curve at 37 GHz we see that indeed there are 
many flares with approximately 4-5 years between them. The 
last flare included in the wavelet analysis peaked in 2002 and an- 
other large flare was observed in 2007, producing a time interval 
of 5 years between the flares. In addition there was a smaller flare 
peaking in 2004 which shows that even though the timescale of 
4.3 years describes the larger events quite well, there is still more 
complex variability present in this source. 

4C 38.41 (1633+382): This is another EGRET-detected 
y-ray bright quasar whic h has been a targ et of many diff^erent 
observing programs (e.g. Katajainen et al. 2000; Jorsta d et al.l 
[2001; Lister & Homan 2005; Bach et al. 2007). In Metsahovi it 
has been monitored for over 20 y ears at 22 and 37 G Hz. It is also 
included in the wavelet study of iKellv et al.l (l2003h which used 
flie UMRAO data at 4.8, 8 and 14.5 GHz. In flieir analysis no 
timescale was found for the source. We found a timescale of 3.4 
years in both frequency bands but the wavelet plot has a complex 
structure and the variability timescale is weaker in the late 1980s 
and early 1990s. The timescale has been present for almost 3/4 
of the flux curve in both frequency bands and repeated 4.8 and 
5.1 times at the 22 and 37 GHz frequencies, respectively. The 
timescale is not seen in DCF analysis of Paper I, probably be- 
cause even though there are some flares with 3.5 years between 
them, it is not a timescale clearly seen to repeat in the flux curve. 

PKS 1749+096: This BLO type object has been monitored 
in Metsahovi for 20 years at 22 GHz and 25 years at 37 GHz. At 
90 GHz the data sets are 14 years long. The source has been a tar- 
get of many Very Long Baseline observations (e.g. iHoman et al.l 
l200ltlWiik et al.l200UlLister & Homanl2005HPiner et al.l2007l) . 
We found a timescale of 2.7 years in all the frequency bands 
but the structure of the wavelet plot is very complex and the 
timescale is increasing in time from 1.7 to 2.7 years. At 90 GHz 
it also repeated only 3.4 times while at 22 GHz 4.4 and at 37 GHz 
5 times. At all frequency bands the timescale has also been 
present for less than 70% of the time. Nevertheless timescales 
of 2.8 years at 22 and 37 GHz and 2.3 years at 90 GHz are also 
seen in the DCF analysis of Paper I. In addition we found flare 
timescales of 0.4 and 0.5 years at 22 and 37 GHz, respectively. 
Both are also seen within 0.2 years in the SF analysis of Paper I. 
The source exhibits rapid variability and has many large flares in 
the flux curve. Indeed many of the flares seem to be within 2.7 
years of each other but there is also other activity present making 
the analysis complex. 

5. Discussion 

Our interest in using wavelets to study the timescales arose from 
the results of Paper I, which showed that many of the sources 
have changed their behaviour during the monitoring time and 
the timescales have changed over the years. A useful property 



of wavelets is that they show also when the timescale has been 
present and how it has changed. Many studies on individual 
sources have claimed the sources to be periodic or quasi-periodic 
based on results from Fourier-based methods. Our results show 
that only a very small number of sources actually show persis- 
tent timescales lasting over the total monitoring time (18 sources 
in which the timescale had lasted for over 90% of the total mon- 
itoring time). A good example is the source 1 156-1-295 shown in 
Fig- [11 The wavelet timescale of 3.4 years is seen to be present 
only in the latter half of the flux curve. In Fig. 1 1 of Paper I 
the DCF and Lomb-Scargle periodogram analyses are shown for 
the same source. The DCF gives a timescale of 3.5 years and 
the Lomb-Scargle periodogram a timescale of 3.3 years which 
are very close to the timescale obtained with wavelets. With 
these methods, however, it was impossible to see that the source 
had changed its behaviour in the mid 1990s and the timescale 
was only present for the latter half of the monitoring period. 
In eight cases the timescale either disappeared or weakened at 
some point during the monitoring period and re-appeared later 
on. These changes cannot be detected by Fourier-based meth- 
ods. Therefore we feel that wavelets should be used more when 
characteristic timescales and quasi-periodicities are studied. 

In many cases the timescale also changed slowly over time. 
In 17 cases we saw a rising trend in the wavelet plot, which 
means a shortening timescale. These were often connected to 
long-lasting flares which showed shorter timescales towards the 
peak of the flare. This could indicate that there are more distur- 
bances in the jet when a flare is growing and more shocks are 
developing in the jet, making the timescale shorter. Also, in six 
cases we detected a declining trend in the wavelet plot, indicat- 
ing an increasing timescale. In these cases the activity is becom- 
ing less frequent and flares do not occur as often as earlier 

In Paper I we also found that the traditional Fourier-based 
methods have many problems when non-sinusoid flux curves are 
studied. In our sample this is the case for all sources. Especially 
with Lomb-Scargle periodogram, we obtain spurious timescales 
and other methods are needed to confirm the timescales. The 
DCF seemed slightly more reliable but it also has no means of 
showing if the timescale is discontinuous or changes in time. We 
compared our results from the wavelet analysis with those ob- 
tained in Paper I. The wavelet timescales are plotted against the 
DCF and periodogram timescales in Fig.[TO| We can see a very 
good correspondence between the diff'erent timescales. We also 
calculated the Spearman rank correlation and between wavelet 
and DCF timescales we obtain a correlation of r = 0.58 {p < 
0.0001). Similarly between wavelet and periodogram timescales 
we obtain r - 0.55 {p < 0.0001). The correlation between the 
wavelet and periodogram timescales is affected by sources that 
have longer periodogram timescales of over 7 years but shorter 
wavelet timescales of 2 to 5 years. In almost all cases the peri- 
odogram analysis also revealed another timescale, similar to the 
wavelet timescale, but in our study we only used the most sig- 
nificant timescale from the periodogram analysis. Only for one 
source (0736+017) did the periodogram analysis give only long 
timescales of over 8 years whereas the wavelet analysis gave 
short 2.7 and 3.4 years timescales. This source, however, has a 
complex structure in the wavelet plot (as seen in Fig.|4]i and mul- 
tiple timescales are seen. 

The average timescales of wavelet and DCF analyses are 
also very similar with the difference being less than 6 months 
at all frequency bands. The averages of the periodogram analy- 
ses are slightly higher The results show that both wavelet and 
DCF analyses give similar results and can be used to determine 
characteristic timescales. The advantage of wavelets is that we 
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Fig. 10. Upper panel: Long-term wavelet timescale against the 
DCF timescale from Paper I. Lower panel: The same wavelet 
timescale against the Lomb-Scargle periodogram timescale from 
Paper I. 



can also take into account how long a certain timescale has been 
present in the flux curve. 

In our sample we had nine sources for which Ke lly et aP 
(I2003h found variability timescales at some of the UMRAO fre- 
quencies. For one of those (1823+568) we were not able to find 
any variability timescale in our analyses. In almost all cases we 
were able to find a similar timescale within 0.5 years of the 
timescale reported at one of the low frequencies. Only in one 
case was the timescale totally different (1641+399) and in one 
case (0212+735 at 37 GHz) there was a difference of 0.7 years 
between the timescales. In three cases it was the timescale we 
refer to as the flare timescale that was similar to the lower fre- 
quency timescale. The similarity between the results from the 
two analyses supports use of the method. 

Even though we found no strict periodicities in the radio be- 
haviour of any of our sources, either in Paper I or in this pa- 
per, we wanted to investigate the possibility of "predicting" up- 
coming active states in these sources. Many of these sources 
have changed their behaviour during the 20 years of observations 
(for example 1156+295 in Fig.[T]), thus such predictions would 
be based purely on their statistical behaviour. This attempt was 
partly motivated by the multifrequency support observations we 
participate in: it would be immensely useful to be able to pri- 
oritise those sources that are likely to enter an active state fairly 



soon. In particular, we were interested to see if we can predict 
which sources would be likely to flare during the early operation 
of th e upcoming garnma-ray satellite GLAST. Our earlier w ork 
(e.g. iTornikoski et alj|2002l: iLahteenmaki & Valtaoiall2003h has 
shown that strong gamma-ray activity usually corresponds to a 
growing flare in the radio regime. Therefore sources that are ac- 
tive in the radio domain are potential sources to be detected by 
GLAST. 

Thus, we wanted to see which sources in our sample would 
be good candidates for strong flaring within the timeframe 2008- 
2009, excluding a set of sources that have almost continuous 
and complex variability at 37 GHz, such as OJ 287 and BL 
Lac. From the wavelet timescales at 37 GHz, combined with the 
timescale analysis from Paper I and with visual examination of 
light curves until the end of 2007, we came up with a list of six 
sources that according to the wavelet analysis and at least one 
other method used in Paper I show semi-periodicity with such 
timescales that a strong radio flare could be anticipated to oc- 
cur within 2008-2009. The sources include one GAL type ob- 
ject (0007+106), two HPQs (0234+285, 1156+295), two LPQs 
(0333+321, 2145+067) and one BLO (1749+096). 

6. Conclusions 

We studied a sample of 80 sources with the continuous wavelet 
transform using data at 22, 37 and 90 GHz. Our aim was to study 
the variability behaviour of the sources and also to better under- 
stand the method and to compare it with Fourier-based methods. 
We found no clear periodicities in the sources. Instead in most of 
the sources the timescales appeared only for a short time in the 
total monitoring period or changed over the years. These kind 
of properties are not revealed with Fourier-based methods and 
therefore wavelets should be used more when quasi-periodicities 
are studied. In our sample only 1/4 had a timescale that had 
been present for at least 4 cycles. In eight of these sources the 
timescale was present in more than one frequency band. This 
shows that the sources are never strictly periodic and conclusions 
about periodicities in the radio frequencies should be drawn with 
caution. 

The average timescale from the wavelet analysis is 4.5 years 
which is the same as was obtained in the DCF analysis of Paper 
I. This timescale should not be affected by the length of the 
monitoring period because we used all the sources in which the 
timescale had repeated at least twice in the analysis and there- 
fore longer timescales also are included in the calculations. The 
weakness of wavelets is that we cannot detect the timescale with 
the same accuracy as with DCF or periodograms and therefore 
both type of methods should be used when periodicities or quasi- 
periodicities are searched for The wavelets should be used to 
verify whether the timescale is persistent or short-lived. If a per- 
sistent timescale is found, the Fourier methods can be used to 
accurately determine the timescale. In our analysis where the 
characteristic timescales are studied, the accuracy of wavelets is 
enough to determine the approximate timescale. 

Acknowledgements. We acknowledge the support of the Academy of Finland 
(project numbers 212656 and 205793). 
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Table 1. Timescales from the wavelet analyses. For each source the long term timescale is given in years and the number of cycles it has repeated. In addition flare timescale and total monitoring 
time are shown. If more than one timescale is given the most significant one is placed first. 
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0007+106 


III ZW 2 


GAL 


19.41 


4.9 


3.5 


0.8 


19.23 


4.9*^ 


3.0 


0.9 


N 


N 


N 


N 


0016+731 




LPQ 


11.87 






0.3 


10.81 








N 


N 


N 


N 


0106+013 


OC 012 


HPQ 


22.27 


6.9 


1.4 


1.7 


24.15 


8.7 


1.6 


2.2 


N 


N 


N 


N 


0109+224 


S2 0109+22 


BLO 


19.33 


8.7 


1.6 




20.66 


2.7'' 


5.0 


1.4 


N 


N 


N 


N 


0133+476 


DA 55 


HPQ 


22.43 


2.2'' 


4.4 




24.17 


2.2'' 


4.4 




14.92 






0.7 


0149+218 




LPQ 


15.89 








16.69 








N 


N 


N 


N 


0202+149 


4C 15.05 


HPQ 


19.93 


4.3*^ 


2.5 


0.9 


20.66 


4.3" 


4.4 


1.1 


N 


N 


N 


N 


0212+735 




HPQ 


15.97 








16.74 


3.4^ 


3.7 




N 


N 


N 


N 


0224+671 




QSO 


13.37 


2.7 


4.6 




14.87 


2.2" 


3.8 




N 


N 


N 


N 


0234+285 


4C 28.07 


HPQ 


16.17 








15.08 


2.4 


4.5 




10.26 


4.3 


1.6 




0235+164 




BLO 


22.44 


5.5" 


3.5 


1.7, 0.3" 


24.11 






0.3 


12.31 


1.9 


2.3 




0248+430 




LPQ 


17.59 








20.73 








N 


N 


N 


N 


0316+413 


3C 84 


GAL 


22.45 








25.46 








14.95 








0333+321 


NRAO 140 


LPQ 


17.28 


4.9 


1.9 


0.8 


25.37 


4.2'' 


3.8 


0.7" 


N 


N 


N 


N 


0336-019 


CTA 026 


HPQ 


15.36 


4.3" 


2.5 


1.7 


16.72 


1.7* 


6.3 




12.28 








0355+508 


NR AO 1 50 


QSO 


22.43 








25.41 








14.95 






0.3,0.7 


0415+379 


3C 111 


GAL 


11.50 


4.3 


2.0 


0.7 


12.46 


4.9 


2.2 




N 


N 


N 


N 


0420-014 


OA 129 


HPQ 


22.32 


5.5 


4.1 




21.11 


5.5" 


2.7 




14.93 


4.3 


2.5 


0.7 


0422+004 


OF 038 


BLO 


19.25 


6.9 


2.4 




19.20 


6.9" 


2.7 


1.7,0.4 


N 


N 


N 


N 


0430+052 


3C 120 


GAL 


22.42 


4.3 


4.4 


1.4,0.5 


24.12 


4.3" 


4.7 


1.4,0.3 


14.95 


2.7 


4.0 


0.7 


0446+112 


PKS 0446+112 


GAL 


15.97 


2.2 


2.5 




16.72 


2.2 


2.0 




N 


N 


N 


N 


0458-020 


PKS 0458-020 


HPQ 


16.16 


6.9 


1.5 




17.08 


5.5" 


3.1 


0.9 


N 


N 


N 


N 


0528+134 


PKS 0528+134 


LPQ 


15.96 








16.92 


3.1 


3.9 




13.47 


2.7 


2.0 




0552+398 


DA 193 


LPQ 


14.04 


6.9 


2.0 




14.96 


8.7 


1.7 


0.9 


13.41 


5.5 


1.5 


0.9 


0642+449 


OH 471 


LPQ 


23.75 






0.9,0.4 


24.10 








N 


N 


N 


N 


0716+714 




BLO 


15.90 


5.5 


2.1 




16.74 


4.3 


1.8 


0.5 


8.82 








0735+178 


PKS 0735+17 


BLO 


23.75 








24.04 








14.94 








0736+017 




HPQ 


21.22 


2.7'-' 


4.0 


0.5 


21.90 


3.4 


4.0 




11.47 


1.4 


4.0 




0754+100 


01 090.4 


BLO 


20.33 


5.5" 


3.1 




25.37 


5.5" 


2.6 


0.9 


N 


N 


N 


N 


0804+499 




HPQ 


15.97 


5.5" 


1.7 


2.2, 0.7" 


16.74 


5.5" 


1.2 


2.2, 0.7" 


N 


N 


N 


N 


0814+425 




BLO 


16.12 


4.3-'^'"' 


2.8 


1.1 


16.89 








N 


N 


N 


N 


0827+243 


0J248 


LPQ 


10.73 




4.2 




11.36 


1.4^-'' 


6.4 




N 


N 


N 


N 


0836+710 


4C 71.07 


LPQ 


16.06 








16.87 


2.7/."' 


4.4 




N 


N 


N 


N 


0851+202 


OJ 287 


BLO 


23.27 


1.4 


10.2 


0.4 


24.82 


1.7 


9.7 


0.3 


14.95 


1.4 


7.2 




0906+430 


3C 216 


HPQ 


20.32 








21.93 








N 


N 


N 


N 


0923+392 


4C 39.25 


LPQ 


23.87 








24.79 








14.91 








0945+408 


4C 40.24 


LPQ 


15.82 








16.73 








N 


N 


N 


N 


0953+254 




LPQ 


15.96 








16.71 








N 


N 


N 


N 


0954+556 


S4 0954+556 


HPQ 


15.82 








16.70 








N 


N 


N 


N 


0954+658 


S4 0954+65 


BLO 


16.14 


5.5"'^'"' 


2.3 


1.4 


21.90 






0.9 


N 


N 


N 


N 


1055+018 


OL093 


HPQ 


22.43 


4.3 


4.4 




24.08 


4.3'' 


4.4 


0.7 


12.74 


2.7 


2.6 




1101+384 


Mrk 421 


BLO 


15.01 








19.13 








N 


N 


N 


N 


1147+245 


B2 1147+24 


BLO 


15.82 








15.79 








N 


N 


N 


N 


1156+295 


4C 29.45 


HPQ 


20.18 


3.4 


1.7 


2.2 


21.92 


3.4 


3.0 


1.7 


12.52 








1219+285 


ON 231 


BLO 


23.89 








24.11 


8.7 


2.5 




N 


N 


N 


N 
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1222+216 


PKS 1222+216 


QSO 


10.55 


2.2 


2.4 


_ 


11.32 


2.2 


2.3 


_ 


N 


N 


N 


N 


1226+023 


3C 273 


LPQ 


24.03 


6.9' 


3.2 


2.2 


25.30 




2.7 


_ 


14.96 


4.3"''' 


3.0 


_ 


1253-055 


3C 279 


HPQ 


24.01 


_ 


_ 


_ 


25.30 


_ 


_ 


_ 


14.93 


_ 


_ 


_ 


1308+326 


AUCVn 


BLO 


22.28 


10.9 


2.0 


_ 


24.00 


9.7 


2.3 


_ 


10.04 


4.3" 


2.3 


_ 


1413+135 




BLO 


15.36 


3.1 


4.5 


_ 


16.13 


3.1 


2.7 


_ 


10.12 


2.7" 


1.8 


_ 


1418+546 


OQ530 


BLO 


21.05 


_ 


_ 


_ 


21.96 


_ 


_ 


_ 


N 


N 


N 


N 


1502+106 


OR 103 


HPQ 


16.11 


6.9"' 


2.0 


_ 


23.47 


6.9"' 


1.7 


1.7 


N 


N 


N 


N 


1510-089 


PKS 1510-089 


HPQ 


19.94 


_ 


_ 


_ 


21.93 


_ 


_ 


_ 


14.05 


3.4" 


2.4 


_ 


1538+149 


4C 14.60 


BLO 


20.33 


2.7c,/,w 


4.0 


_ 


21.93 








N 


N 


N 


N 


1606+106 


4C 10.45 


LPQ 


11.28 


5.5 


1.5 


_ 


12.10 


4.9 


1.9 


_ 


N 


N 


N 


N 


1611+343 


DA 406 


LPQ 


15.96 


4.3,10.9 


3.1 


_ 


16.83 


4.3,10.9 


3.1 


_ 


N 


N 


N 


N 


1633+382 


4C 38.41 


LPQ 


22.43 


3.4^ 


4.8 


_ 


24.08 


3.4" 


5.1 


_ 


N 


N 


N 


N 


1637+574 


OS 562 


LPQ 


20.17 


3.4'' 


4.0 


_ 


21.94 


3.9'^ 


2.8 


_ 


N 


N 


N 


N 


1641+399 


3C 345 


HPQ 


24.00 


8.7'- 


2.5 


- 


24.80 


8.7'- 


2.7 


- 


14.93 


- 


- 


- 


1652+398 


Mrk501 


BLO 


16.11 


- 


- 


- 


16.91 


- 


- 


- 


N 


N 


N 


N 


1725+044 


PKS 1725+044 


QSO 


12.70 


- 


- 


- 


13.62 


- 


- 


- 


N 


N 


N 


N 


1739+522 


S4 1739+52 


HPQ 


15.84 


5.5 


2.4 


1.1 


16.87 


4.3'^''' 


2.5 


1.1 


N 


N 


N 


N 


1741-038 


PKS 1741-038 


HPQ 


16.03 


5.5'- 


2.8 


- 


16.99 


5.5'- 


2.9 


0.5 


13.35 


6.1 


2.0 


- 


1749+096 


PKS 1749+096 


BLO 


19.86 




4.4 


0.4 


24.53 




5.0 


0.5 


14.10 




3.4 


_ 


1803+784 


S5 1803+784 


BLO 


15.91 


3.4^-"' 


2.9 


_ 


16.75 


3.4^'"' 


2.9 


_ 


9.10 


_ 


_ 


_ 


1807+698 


3C 371.0 


BLO 


20.17 


_ 


_ 


_ 


21.95 


_ 


_ 


_ 


N 


N 


N 


N 


1823+568 


4C 56.27 


BLO 


15.56 


_ 


_ 


_ 


16.70 


_ 


_ 


_ 


9.04 


_ 


_ 


_ 


1928+738 


4C 73.18 


LPQ 


16.06 


5.5"' 


1.9 


_ 


16.85 


_ 


_ 


0.6^ 


N 


N 


N 


N 


2005+403 




QSO 


22.29 








24.03 








N 


N 


N 


N 


2007+776 


S5 2007+77 


BLO 


12.38 


3.1"' 


2.7 


_ 


16.83 


_ 


_ 


_ 


N 


N 


N 


N 


2021+614 


OW637 


LPQ 


16.80 


. 


. 




21.93 








N 


N 


N 


N 


2134+004 


OX 057 


LPQ 


22.27 


6.9 


3.2 




25.37 








N 


N 


N 


N 


2136+141 




LPQ 


15.51 


8.7 


1.5 




18.20 


6.9 


1.4 




N 


N 


N 


N 


2145+067 




LPQ 


18.37 


6.1 


2.7 


- 


19.16 


5.5 


3.0 


- 


14.93 


- 


- 


- 


2200+420 


BLLac 


BLO 


24.01 


8.7'- 


2.5 




25.44 


l.T 


3.1 




14.96 






0.9 


2201+315 


4C 31.63 


LPQ 


19.25 


4.3 


2.5 




22.54 


A3 


4.4 


0.5 


N 


N 


N 


N 


2223-052 


3C446 


BLO 


19.35 


9.7 


1.9 




19.22 


9.7 


1.8 




14.94 


8.7 


1.5 




2230+114 


CTA 102 


HPQ 


19.25 


8.7 


1.9 


0.8 


19.21 


6.9 


2.7 


0.9 


14.50 


1.1 


1.6 




2234+282 




HPQ 


15.97 








15.83 








N 


N 


N 


N 


2251+158 


3C 454.3 


HPQ 


24.01 


6.P 


3.6 




24.48 


6.P 


3.6 


1.4 


14.95 









" = complex structure in the wavelet plot. 

= weak trmescale. 
f = faint source. 
'' = rising trend in the timescale. 

= declining trend in the timescale. 
- = timescale not determined. 
N = not enough data for wavelet analysis. 
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1985 1990 1995 2000 
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Fig. 2. Wavelet transform of the source DA 55 (0133+476) at 37 GHz. 
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Fig. 3. Wavelet transform of the source 3C 120 (0430+052) at 37 GHz. 
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Fig. 5. Wavelet transform of the source OJ 248 (0827+243) at 37 GHz. 
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Fig. 6. Wavelet transform of the source OJ 287 (085 1 +202) at 37 GHz. 
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Fig. 7. Wavelet transform of the source OL 093 (1055+018) at 37 GHz. 
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Fig. 8. Wavelet transform of the source 4C 38.41 (1633+382) at 37 GHz. 




1985 1990 1995 

year 

Fig, 9. Wavelet transform of the source PKS 1749+096 at 37 GHz. 
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